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SUMMARY 


A theoretical analysis of the effect of transverse shear and 
rotary inertia on the natural frequencies of a uniform "beam is 
presented. Frequency equations are derived for the cases of the canti- 
lever heam, the sy mme trically vibrating free-free beam, and the anti- 
symmetrically vibrating free-free beam. Numerical results are pre- 
sented in the form of curves giving the frequencies of the first three 
modes of the cantilever beam and the first six modes, three symmetrical 
and three anti symmetrical, of the free-free beam. 


INTRODUCTION 


In the dynamic analysis of aircraft, the determination of the 
natural frequency is of basic importance. A number of methods of 
vibrational analysis have been presented in the past. Tests cm actual 
structures, however, have often shown discrepancies between the calcu- 
lated and the observed values of the natural frequency. 

One of the possible explanations for these discrepancies is that 
many of the previous calculations were based an methods in which only 
the elementary engineering theory of beam bending was used and no 
secondary effects were included. Among the secondary effects are shear 
lag, deformation of the webs due to transverse shear, and rotary 
inertia. The effect of shear lag and shear deformation of the web is 
to increase the flexibility of the beam because of the additional 
deflection that is Introduced. The effect of rotary inertia is to 
increase the dynamic loading on the beam because of the additional 
inertia loading due to the rotational acceleration of the differential 
elements of the beam. 

The effect of shear lag on the bending vibrations of box beams was 
discussed in reference 1 and is not considered herein. The effects of 
transverse shear and rotary inertia were discussed by Timoshenko 
(reference 2) but only to the extent of presenting the differential 
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equation. In the present paper a general solution of this differential 
equation is given together with specific solutions for the cantilever 
and free-free beams. (See appendix.) Charts giving the frequencies of 
the first three modes of the cantilever and the first six modes (three 
symmetrical and three anti symmetrical) of the free-free beam are 
presented. 
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SYMBOLS 

1 shear stiffness, pounds 
shear area, square inches 

effective total cross-sectional area, square Inches 


Ci, c 2 , c 3 , C 4 

El 

E 

G 
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L 

M 
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constants of integration 

flexural stiffness, pound- inches 2 
modulus of elasticity, psi 
shear modulus, psi 
moment of inertia, inches^ 

length of cantilever beam and half-length of free- 
free beam, inches 

moment, pound- inches 

reversed effective moment due to rotational accelera- 
tion, pound-inches 

shear, pounds 

acceleration due to gravity, inches per second per 
second 

frequency coefficient. 



frequency coefficient where shear and rotary inertia 
are neglected 
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coefficient of shear rigidity 



coefficient of rotary inertia ( i \|-i- 



p p 

mass of "beam per unit length, seconds per inch 
distributed loading, pounds per inch 
time, seconds 

weight of beam per unit length, pounds per inch 

distance along span measured from root or center 
line, inches 

deflection of beam, inches 
constants defined with equation (5) 
shear strain 

rotation of cross section, radians 

natural frequency of beam, radians per second 

natural frequency of beam excluding secondary effects, 
radians per second 

nondimensional crordinate 

density of beam, pounds per cubic inch 


RESULTS AND DISCUSSION 


Theoretical .- The natural frequencies for the cantilever beams and 
the free-free beams are defined by the frequency equations derived in 
the appendix. These equations are in terms of the parameters kg, kg, 

kgj. In order to solve the equations, a trial-and-error process 

and kc- there are 


is used. For every combination of values of kg 
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an Infinite number of values of kg that satisfy the frequency 
equations. The smallest value of kg is associated with the first 
mode, the next larger value with the second mode, and so forth. 

Numerical .- The frequency equations were used to calculate natural 
frequencies for the first three modes of the cantilever and the first 
six modes (three symmetrical and three antisymmetrical) of the free- 
free team. The results of these calculations are given in c har t form 
in figures 1 to 3* In these charts the ratio of the natural fre- 
quency cd to the natural frequency coq obtained ty neglecting second- 
ary effects is plotted as a function of the shear- stiffness param- 
eter kg and the rotary- inertia parameter k RI . In order to obtain cd 

from the ratio cd/cd d , cd q must first be calculated from the formula 


where k B 0 is a constant for a particular mode and is given on the 

chart for that mode in figures 1 to 3. The values of cd/cd q are given 
only for values of kg from 0 to 0 . 2 k and for values of k RI in 

increments of 0 .05 from 0 to 0 .20 inasmuch as. values above these are out 
of the range of probable design. The frequency equations derived in 
the appendix, however, can be used if values of cd/cd q are needed for 
values of kg or k R j greater than those shown in the figures. 

Inspection of the curves shown in figures, 1 to 3 shows that in 
extreme cases for which the coefficients k s and kpj are large, 

reductions in frequency as hirdi as 50 percent can be had in the first 
mode; whereas reductions as high as 80 percent can be had in the third 
mode. , 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Ya., May 1 6, I9I4-9 
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APPENDIX 

DERIVATIONS AND SOLUTIONS OF FREQUENCY EQUATIONS 


In the solutions for the natural tending frequencies of teams 
tased on the elementary engineering theory, the team deflection is 
considered as a function only of the manner of support nnfl the flexural 
stiffness El and, also, dynamic loading is considered as a function 
only of the translatory acceleration of the particles. The deflection, 
however, is also influenced ty the shear stiffness AgG of the team 

and the dynamic loading is influenced ty the rotational acceleration of 
the cross section of the team. 

The effect of the shear deformation is introduced into the solu- 
tion in the equation for the slope of the deflection curve. In the 
engineering theory of tending, the slope of the deflection curve 

dx 

is assumed to depend only on the rotation of the cross sections of the 
team associated with flexure 0 . The slope, however, depends also an 
the shear strain 7 which when included gives 


d£ 

dx 


0+7 




where M is the moment and V, the shear. The coordinate syBtem is 
shown in figure 4. 

The effect of the rotary inertia is introduced into the solution 
in the summation of moments of the differential element shown in 
figure 5 • The elements of the team perform not only a translatory 
motion tut also a rotation. The rotation 0 is due to tending alone 
since the shear just superimposes a sliding of adjacent cross sections 
with respect to one another. This rotation varies with time. The 
angular velocity and the acceleration of each element are given 

ty ^ and respectively. The reversed effective moment of the 

dt2 

forces due to this angular acceleration on a length dx is expressed as 
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where p is the density of the "beam and g is thet acceleration due to 
gravity. The moment is taken as positive in the clockwise direction as 
shown in figure 5 • 

The inclusion of both of these effects in the equilibrium equa- 
tions of a uniform beam in free harmonic vibration results in the 
following differential equation: 


(” 2 - ■ 0 w 


At effective total cross-sectional area (w/p) 

w weight of the beam per unit length 

m mass of the beam per unit length 

to natural frequency 

The terms not containing Ag and Ap are the terms that make up the 

standard differential equation when the effects of shear stiffness and 
rotary inertia are neglected. The middle term is composed entirely of 
secondary effects . The last term contains an interaction effect of the 
two secondary effectB. 

Division of the differential equation by nu>2 and the substitution 
of the nondimensional coordinate I = — change equation (1) to 


where 


El 




mcp2 

AsG 



effective shear-carrying area 


1 d' 


k B 2 di 


<£y 


( k s 2 


■ ( 1 ‘ " 0 


( 2 ) 


where 



(3a) 
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(31>) 


and 




(3c) 


The quantity can he thought of as a frequency constant, as can he 

seen hy transposing equation ( 3 a) to 

W 


For a given uniform "beam the quantity ip^r- is constant; therefore, 

j|mL 4 

the natur al frequency of* the teem is a function of the coefficient k B . 
The solution of the differential equation can he written in the • 

form 

y = Ci cosh k B a£ + C 2 sinh k B a4 + C 3 cos k B P| + C 4 sin kgp| (5) 


where 


a = 


I 


-(kg 2 + k EI 2 ) + ^(kg 2 - k RI 2 ) + 




P = 


(kg 2 + k EI 2 ) + ^(kg 2 - %X 2 ) 2 + 


and C^, C 2 , C 3 , and are constants of integration. 
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Boundary conditions The "boundary conditions for the vibrating- 
beam problems considered in this paper will he specified at the root 
and at the tip. 

For the conditions at the root, the equations are as follows: 

For the cantilever team the conditions of zero deflection and 
slope equal to shear strain are expressed by 

M|=o = 0 ^ 


ilk 

I«[d|_ 


1=0 




k S 2fc B 2 


1 - k S 2k RlV 


L^Z + 

L^b 2 d ^ 3 



( 7 ) 


For the symmetrically vibrating free-free beam, the conditions of zero 
shear and, therefore, zero slope are expressed by 


d£| = 0 (8) 

dy |=o 



( 9 ) 


For the antisymmetrically vibrating beam the conditions of zero deflec- 
tion and zero moment are expressed by 



d 2 y 

_dt 2 _ 


= 0 


(10) 


1=0 


( 11 ) 
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For the conditions at the tip for all beams, the conditions of zero 
moment and zero shear are expressed by the following equations: 



The complete deflection curve of a "beam is obtained by substi- 
tuting equation (5) in the proper boundary conditions. This substi- 
tution results in a set of four homogeneous linear equations 
in C 2 , C^, and . In order that solutions other than zero 

exist (that is, that vibration can occur), the determinant of the 

coefficients of the C's must be equal to zero. 

Frequency equations .- The equations obtained by setting the 
determinant equal to zero give the relationship between kg, k RI , 

and kg required to determine the natural frequency which is contained 

in k-g. These equations for each beam treated in this paper are as 
follows : 

For the cantilever beam 

„ k B( k S 2 + k RI 2 ) 

2 ‘ 1 o- = 5 == g 8ln kBp Sinil k B a 

V 1 - *S k RI k B 

k “ k pj 2 ) I cos k-gP cosh kga = 0 ( 14 ) 
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For the symmetrically vibrating free-free beam 


P^P 2 - kg 2 )tanh kg<x + a (a 2 + k g 2 ) tan kgP = 0 ( 15 ) 


For the antisymmetrically vibrating free-free beam 


a(a 2 + k g 2 )tanh k B a - p(p 2 - kg 2 ) tan k B p = 0 ( 16 ) 
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Figure 5. — Forces on a differential element. 




